On the multifractal structure of fully developed turbulence. 
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A new way to introduce stochasticity into the Navier-Stokes equation is proposed. On this basis, 
the power-law dependence of velocity and vorticity correlators and their multiscaling properties are 
derived from the Navier-Stokes equation. 
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INTRODUCTION 

Understanding of statistical properties of a turbulent 
flow is a classical problem of hydrodynamics. Intermit- 
tent behavior of velocity scaling exponents demonstrated 
in experiments and numerical simulations [1-3] is one of 
its important aspects. The most successful, and conven- 
tional nowadays, way to interpret this intermittency is 
the Multifractal (MF) approach introduced in [4]. This 
is a generalization of Kolmogorov's K41 theory, and it 
allows to express all the observed intermittent charac- 
teristics by means of one function D{h). In particular, if 
D{h) is derived from observed velocity scaling exponents, 
then all the other values can be calculated [5]. 

However, the phenomenological presentation of the MF 
model assumes the existence of singularities (actually, the 
paper [4] was titled 'On the singularity structure of fully 
developed turbulence'). The question if singularities can 
appear after a finite time has been widely discussed and 
remains still open [6-8] . But in the absence of even one 
example of a finite-time singularity, the existence of a 
whole 'spectrum' of singularities is rather doubtful. 

To make the MF theory independent of this as- 
sumption, its probabilistic reformulation was introduced 
[6, 9, 10]. However, this formulation does not allow to un- 
derstand what happens in a turbulent flow; the structures 
(or the solutions to the Navier-Stokes equation) that are 
responsible for the observed intermittent properties, re- 
main unknown. 

On the other hand, many experiments and numerical 
simulations have shown the presence of 'coherent struc- 
tures' - vortex filaments and 'pancakes' - in a turbulent 
flow (see, e.g., [11]). In [12] it was shown that these struc- 
tures make a fundamental contribution to the observed 
Kolmogorov's two-thirds law, and they contain the most 
part of the whole enstrophy of a flow. The origin of these 
coherent structures is still obscure (in [6] the tendency of 
a flow to produce these structures was characterized as 
'mysterious'). Understanding of these filaments forma- 
tion seems to be very important as it may also help to 
understand the cause of intermittency and multifractal- 
ity. 

In our previous papers we have developed a model of 
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vortex filaments (VF). We showed that it did not con- 
tradict to the MF model [13]. In [14, 15] we combined 
the VF theory with the MF approach; this allowed to 
calculate both longitudinal and transverse velocity scal- 
ing exponents without adjusting parameters, the result 
agreed well with numerical data. Although this model 
was based on the NSE, it still contained an additional 
supposition: it assumed that the nonlinear part of the 
pressure hessian was orthogonal to the local vorticity in- 
side the vortex filaments. This assumption was confirmed 
by reasonable considerations but it was not proved thor- 
oughly. 

In this paper we propose a new accurate approach that 
may help to derive the multifractality and the existence 
of vortex filaments based on the Navier-Stokes equation 
(NSE). In our approach, the power spectrum appears as 
a result of stretching of a vortex filament. Scaling ex- 
ponents of different orders are produced by different fil- 
aments contributing mostly to these orders. No finite- 
time singularities are needed, singularity is reached (in 
non- viscous limit) after infinite time. The power-law de- 
pendence of velocity in the vicinity of the forming sin- 
gularity (not the cascade of decaying vortices) produces 
the long power-law tail of the Fourier spectrum. 

A new way to introduce randomness into the NSE 
or Euler equation allows us to consider the evolution of 
small-scale velocity perturbations in the 'external' field 
produced by random large-scale velocity (Section 2). We 
find a long-time asymptotics of general solution for small- 
scale fiuctuations (Section 3). It appears that this asymp- 
totics depends on some combination of large-scale val- 
ues, which is random, but tends to a constant as time 
increases. Analyzing the solution, we observe an effect 
similar to that described in [6] on the basis of numerical 
simulations: in some special reference frame, velocities 
depend to leading order on only one coordinate, thus the 
fiow becomes one-dimensional; this causes a depletion of 
nonlinearity. The evolution along this one coordinate 
corresponds to stretching of the rotating vortex filament. 

In Section 4 we thus introduce a simplified model, in 
which the large-scale field is fixed and the special ref- 
erence frame coincides with the laboratory frame. This 
particular case allows to understand better the details of 
solution, to see what happens to the spectrum, and to 
estimate the effect of viscosity. 

In Section 5 we discuss the first-order correction to 
the solution; not only averages of the combined large- 
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scale value mentioned above, but also fluctuations around 
them are now taken into account. We see that fluctua- 
tions of the large-scale flow result in multifractality of 
small-scale correlators. 

In the last section we discuss the results of the paper. 

EQUATION FOR SMALL-SCALE VELOCITY 
FLUCTUATIONS 

A turbulent flow can, in principle, be completely de- 
scribed by the dynamical NSE. The pulsations appear 
as a result of instability of the flow. But solving the 
complete problem with account of instabilities seems to 
be impossible. So, to derive statistical properties of the 
flow, one has to introduce randomness into the equation. 
This is usually done by adding a large-scale random ex- 
ternal force into the right-hand side of the NSE: 

_ + (vV)v = -VP + F{r,t)+iyAv, (1) 
V-v = 

The probability distribution of the force is usually sup- 
posed to be Gaussian. One assumes that the resulting 
correlation properties inside the inertial range do not de- 
pend on the properties of the exciting force. To satisfy 
this requirement, F must not include small-scale pulsa- 
tions: not only its correlator must decay at scales of the 
order of largest eddies' turnover scales L, but also the re- 
alization of F must consist of large-scale harmonics only. 

But the external volume-acting forces might exist in 
the flow and might not. On the other hand, this way does 
not allow to separate; large-scale and small-s(;alc velocity 
fluctuations; introducing the large-scale force does not 
make any simplification. 

We propose a new approach to the problem: the large- 
scale velocity perturbations are considered as given ran- 
dom process. The stochastic properties of small-scale 
fluctuations can then be derived based on the properties 
of the large-scale fluctuations. 

Formal introduction of randomness 

To develop this idea, introduce a random field U. 
To make it large-scale in the sense discussed above, we 
smoothen U by means of a space average over surround- 
ing volume of the order of L, e.g.: 

Ui(r, t) = ^JUi{v + p, t)e-f'/'^'dp (2) 

and we require V • U = 0. 
Now we define the large-scale force F according to: 

-^ + (UV)U = -V7r + F(r,f)+i/AU, (3) 
V-F = 



It is evident that F is large-scale and satisfies the above 
condition. (To avoid difficulties with the time derivative, 
one should consider weak solutions wherever it is needed; 
as wc will sec later on, only U and its integrals, not 
derivatives, contribute to the result.) 

We then substitute this F to the right-hand side of (1), 
and seek the solution in the form 

v(r,i)=U-Fu, P = p + 'JT 

Then for u we get the equation: 

§^Ui + (UV)Mi + (uV)[/i + (uV)Mi = -Vip + lyAui , 
ViUi = (4) 

This is the stochastic version of the NSE with stochas- 
ticity introduced by means of the large-scale random ve- 
locity field U (instead of random force). 

We note that F is not presented in this resulting equa- 
tion. So, in this approach it is not important if the ex- 
citing force docs exist or not. In this sense, introducing 
probability by means of the large-scale velocity field U is 
more general and includes the NSE with external forces 
particular case. 

As it was with large-scale forces, we assume that any 
particular choice of U does not aff'ect significantly the 
statistics of the flow at small scales. Indeed, as wc will 
see below, the choice of the statistical properties of U, 
and in particular its gaussianity, is not important for the 
result. 

Physical reformulation 

In our approach, U can be treated as some (large-scale) 

part of velocity field that is chosen to be an independent 
random process. One could as well define the whole large- 
scale component of v to play this role. To split a velocity 
perturbation into large-scale and small-scale components, 
one can take a space average of v over surrounding vol- 
ume of the order of L: 

Vi{r, t) = ^ Jvi{^ + P, t)e-p""'"dp , V • y = (5) 

Pressure and forces can be smoothed in the same way. 
Applying the 'smoothing' procedure (5) to the NSE and 
neglecting the contribution of the 'smooth' part of the 
non- linear small-scale term at large scales,^ one gets the 



[1] Actually, if v = V + u, the contribution of small scales to the m- 
erage is J(u ■ V)«i(r -|- p, t)e~'' dp ^ J Ui\i- dS, where 
the integral is taken over the sphere of radius L. The expectation 
value of this expression is zero; thus, one can estimate the term 

by -j-^ ij D(J UiU ■ dS) oc -j^. Here wc neglect the contribution 
of fluctuations at intermediate scales (see Discussion). 
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equation describing the evolution of the large-scale terms. 
The difference between this equation and the NSE gives 
Eq. (4) with U = V. In accordance with the above 
idea, one now treats V as independent random function of 
time. The equation then becomes a stochastic equation. 

Thus, (4) can as well be treated as the equation for the 
small-scale component of velocity u, while the large-scale 
component of the same velocity field acts as a source of 
randomness. 



where gin{t) and qva{t) satisfy the equations: 



9ia ~t~ -^ijOjoi ? 

Substituting to (7) we get: 

' dw,. 



qij{0)=5ij 



dw,, 



dp 



(9) 



The small-scale limit 

We now simplify (4) to analyze it analitically. The 
smoothed function U can be expanded in a Taylor series 
for r ^ i, so 



Ui (r, t) = C/i (0, t) + {t)rj + . . . , Ai, 







(6) 



The omitted terms are smaller by a factor ^ r/L. Wc 
hereafter restrict our consideration to the first two terms 
since in what follows we will be interested in small scales. 
This corresponds to taking the limit L — )■ oo with large- 
scale eddy turnover time T remaining constant. Note 
also that the drift component U(0, t) can easily be taken 
zero by choosing the appropriate reference frame. 
Then from (4) we obtain: 



d 

-g^Ui + {AkjrjVk)ui + AikUk + {uV)ui 



ViUi = 



(7) 



This is the main equation of the paper. In the limit 
r <^ L, it is an exact consequence of the Navier-Stokes 
equation. The large-scale velocity gradients Aij play the 
role of external forces, and, instead of the forces, they are 
used to introduce stochasticity into the equation. Know- 
ing their statistical properties, one can now analyze the 
equation to find statistical properties at small scales. 

The experiments show that large-scale fluctuations are 
close to Gaussian [6]. Since cutting the series (6) can 
be treated as taking the limit L — >■ oo, the Gaussian 
approximation is the more confident. However, for our 
purposes the gaussianity is not necessary of A^j . 

Our task is now to study the asymptotic properties of 
the stochastic equation (7) after long time. 



ASYMPTOTIC ANALYSIS OF (7): INVISCID 
LIMIT 



Wc now consider the Euler analog to (7), i.e. v = 
(the contribution of viscosity will be discussed later). 
First, we change the variables r, u to X, w according 

to 



dwfj. 



= 



In this paper we restrict ourselves by the consideration of 
symmetric Aij-. we will discuss the inner parts of vortex 
filaments where vorticities are very high, so we expect 
that small 'external' large-scale vorticity (which is equal 
to the asymmetric part of Aij) does not play a crucial 
role. Prom A = A^ it follows 

9ij — Qji 



The equation then becomes 

dt + 1-^9jaWa g-^^ 



dp 

"dx: 



(10) 



dX, 



Since A^j is a random process, the matrices gij and qij 
are also random. Note that only the combination qjigik 
is presented in (10). To analyze the solutions of (10) at 
t — >■ 00, we need to know the asymptotic behavior of this 
value. 



Asymptotic behavior of q,g and eisymptotic solution 
of (10) 

To examine the solution of (9), we proceed to a discrete 
approximation: consider a discrete sequence of moments 
separated by At and let Aij{t) = (A„)jj be constant 
inside each small (n-th) interval. Then, for each Af, the 
solution to the Eq. (9) is described by an exponent, and 
we get 



qn = qn-ie 



Hence, 



We now consider the Iwasawa decomposition of the ma- 
trix qj^: 



Ui{r, t) = giij.{t)w^{X^, t), X^ = q^a{t)ra (8) 



q = z{q)d{q)s{q) 



(12) 
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where z is an upper triangular matrix with diagonal el- 
ements equal to 1, d is a diagonal matrix with positive 
eigenvalues, s is an orthogonal matrix. 

The matrix qjsi is a multiplication of N random real 
unimodular matrices with the same distribution. The 
asymptotic behavior of this object has been studied care- 
fully, and a number of important results has been ob- 
tained. (For short summation of them, see [16].) In par- 
ticular, the following Theorems have been proved under 
reasonable conditions: ^ 

1. [17] with probability 1, there exists the limit 

lim iln(ii(gAr) = Aj, Aj is not random, and 

Ai < A2 < A3; 

2. [18, 19] the distribution of $i = (inrf.(9»)_-A.iV) 
asymptotically close to a Gaussian distribution and 
(weakly) converges to it as n — )• 00; 

3. [20] with probability 1, 2;(9jv) converges as A'' — )• 00, 

2:(gjv) 

4. [21] the values ^i{qN) and ^(gjv) are asymptotically 
independent. 

For our purposes, these results can be written shortly 
3 



as 



z{Qn) Zoo , 

N N N 

, , , AiiV-KECi(") ^2N+j:i2in) \3N+j:i3{n) 

d{qN) = diag(e " ,e " ,e 



Ai < A2 < A3 



(13) 



(The ordering is due to the triangular matrix.) Note 
that the matrix z^o and the coefficients Aj are constants 
determined by statistical properties of the random pro- 
cess (An)ij. To the contrary, s^qn) changes quickly as a 
function of N, and depends strongly on the realization of 

In the case A = , fortunately, the rotating term 
vanishes in the combination qg: 

i<19)N = {qq^)N ^ Zood{qN)z'^ , 

d{qN) = e^^^^ • diag(0, 0, 1) + 0(e^^^) 
Neglecting the terms growing slower than e^^^'^ , we get 
{qg)N = Ce'^'" 



[2] Note that the Theorems assume neither symmetry nor gaussian- 
ity of A. 

[3] The transition to the hmit N 00 corresponds to t — > 00 with St 

remaining constant. To return back from discrete to continuous 
description, one lias to take the limit At — s- afterwards. This 
would result in replacement A'^ by t and rcnormalization of Aj. 
The exponents in (13) would take the form Ait -|- / ^idt. 



where, with accuracy O (^e^^^~^''^^) , C is a constant sym- 
metric matrix C = Zoodiag{0,0, 1)2;^. 

We now introduce a new vector variable U = Cw in- 
stead of w, and wc return to the continuous description 
qg{t). Then from (10) wc get 



dV 
'dt 



o2A3t 



u 



_d_ 



U = -C 



dP 



dU 

ax 



= 



From incompressibility it follows Ai -|- A2 -t- A3 = 0, hence 
A3 > 0. Thus, the asymptotic {t 00) solution to the 
first order by e~^^^* takes the form: 



U 



_d_ 
5X 



5X 



0, P 



o2A3t 



n 

(14) 

This equation is the asymptote of Eq.(lO) at large values 
of t. We note that the constant matrix dj is the only 
remnant of the random process Aij{t) in this equation. 
This is of course due to the chosen variables (X. U): the 
randomness remains in rotation of the corresponding ref- 
erence frame. 

The matrix C is symmetric and can be reduced to di- 
agonal by some (time-independent) twist of the reference 
frame. Thus, the solutions of (14) correspond to some 
stationary hydrodynamical configurations. 

So, we see that the relations (8),(9),(12)-(14) set up 
a correspondence between general stochastic solutions of 
Eq. (7) at very large values oi t {t ^ 00) and stationary 
solutions of the Euler equation. 



Analysis of the solution 

To understand the properties of the solution (14), we 
have to rewrite it back in laboratory coordinates (r, u). 
With account of (8), (9), (12) we have 

u = gw = gC-^U{X) = q^C-^U{X) = s^dz^C-'^U{X) , 



X = gr = zdsr 

To separate the stochastic rotational part of the solution, 
we make one more change of variables: 



sr 



u = su 



(15) 



This frame rotates randomly, since the matrix s is a 
stochastic function of time (as opposed to z and d, which 
tend to a constant or change steadily at large t). Also 
for any U(X) determined by (14), we define a new vector 
function 



Then 



V(y) = z^C-'U{zy) 



u' = dV{dr') 
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or, in more detailed writing, 

u'^ = e^*Vi(e^iV'i,e^=*r^,e^^V^) 

(no summation is assumed). 

We see that in the rotating coordinates r', the asymp- 
totic solution is not random. 

As t ^ 00, the third component u'^ dominates, and the 
solution stretches exponentially with different cocfScients 
along different axes. Hence, to the leading order it is 
enough to account the dependence of u on only one {r'^) 
variable.'* 

We now take the curl to find vorticity: 

Since ^ = 0, we have cjk « e~'^'=*. Hence, vorticity is 
directed mainly along the r[ axis: 

u;':^Lu[= e-^^V (16) 

We note that, since a;' = suj, the absolute values of vor- 
ticities are equal in the two frames, so w = w'. 

Thus, vorticity (and velocity) is transported from 
boundaries to the center, and simultaneously it grows 
exponentially. To keep the whole system stable, we have 
to demand that at some point, e.g. r'^ = L, vorticity (or 
velocity) is nearly constant: 

w(t,L)~l (17) 

In fact, the requirement is much weaker; it would be 
enough for uj(t, L) not to grow exponentially. This is true 
for a point of general position.^ However, here we ask w 
to be nearly constant at the boundary for simplicity. We 
shall discuss the subject in the next Section. 

With account of the boundary condition, we have 
/(e^^* L) e^i* for any t'. Choosing t' in such a way 
that e^^^r'^ = e"^^* L, we rewrite the solution (16) in the 
form: 

^{tys)^[j) (18) 

The conditions of the Theorems on page 4 require t to be 
large enough, so the equation (13) is valid for some t > to- 



[4] Prom the condition V • u = it follows 

— • u' = Ve2^**^(e^i*ri,e^^V^,e^3V') = 
dr' ^ dyi ^ 

[5] Indeed, if the measure of the points where vorticity grows expo- 
nentially with time is ^ > 0, then velocity also grows exponen- 
tially in the region. Then, to satisfy the condition that energy 
docs not grow exponentially on average, // must decrease expo- 
nentially itself. 



Then (18) is valid for all t > t*{r'^) = ^ In(ryi) + to, 
or r'^ > Le^^^^°~*^\ At smaller r'^, the influence of the 
boundary has not yet reached the region, and w is deter- 
mined by the initial condition. So, (18) does not mean 
a real finite-time singularity: it is naturally 'smoothed' 
near the center, the radius of the smoothing part steadily 
decreasing. On the other hand, (18) gives a power law 
for velocity structure functions. Scaling (not yet multi- 
scaling) properties are derived from the stochastic Euler 
equation. 

We stress that (18) is a long-time local approximation 
to any general solution (14) of non-viscous Eq. (7) with 
symmetric Aij in the regions of very high vorticity. In 
the next Section we discuss a simplified model that helps 
to understand the solution better and to reveal the role 
of viscosity. Then we proceed to the discussion on mul- 
tifractality. 



SIMPLE MODEL: DETAILS OF SOLUTION AND 
ACCOUNT OF VISCOSITY 

In the previous section, the long-time solution of (7) 
was found to be nearly one- dimensional in the rotat- 
ing frame (15), and its behavior appeared to be rather 
deterministic. The randomness is contained mostly in 
the rotation matrix s. Here we propose a simple one- 
dimensional deterministic model equation that has simi- 
lar solutions and helps to understand the details of their 
behavior. Later on, we use this model to generalize the 
results for finite viscosity and for multifractal description. 

The idea is to 'straighten' the random flow, excluding 
the matrix s and thus avoiding the need of additional 
rotation (15) of the frame. So we fix the random matrix 
Aij and restrict ourselves by small-scale velocity field de- 
pending on only one variable. 

Thus, consider the velocity field: ® 

Vx=a{t)x, Vy = b{t)y + u{x,t) , = c{t)z 

The parameters a, b, c correspond to the large-scale ma- 
trix Aij. From incompressibility it follows 

a+b+c=0 

The Euler equation then takes the form 

^(« + «')=-|f 

^ + a{t)x'-^ + b{t)u{x, t) + y(b + b^)=-% 
-(c + c^)=-if 



[6] More general consideration with u = u(a;, y, z) gives the same 
results. 
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The pressure derivatives must depend linearly on x,y, z, 
respectively. Hence, pressure must take the form: 

The values pi,p2,P3 are determined by the evolution of 
a, b, c, respectively; the equations for them are equivalent 
to (3). The part of the second equation that does not 
depend on y is: 



dt 



dx 



This equation describes the evolution of the small-scale 
component and is analogous to (7). 

In the region of interest velocities are small, while vor- 
ticities are very high. So, in what follows we will discuss 
vorticity instead of velocity (although very similar rela- 
tions can be written for velocity). Since uj{t,x) — 
the corresponding equation is 



- + a(t).--cWa. = 



(20) 



We will hereafter analyze the solutions to this equation 
in the range x € [0, 1] for t > 0. For simplicity, let a, h 
be constants (although the solution can be written for 
arbitrary functions a{t),b{t)). Let also 

a<0, b>0, c=-{a + b)>b 

In addition, we set the boundary condition 

w(t,l) = l (21) 

(This is a strengthened variant of (17).) To provide this 
boundary condition, the initial condition w(0, x) — ujo{x) 
must satisfy 

c.o(l) = l, a^(l)-c = 

It is easy to check that all solutions of (20) obey the 
relation: 

a;(i, x) = e"(*-*')a; (t' , xe-'^f*"*'') (22) 

For any x > e"^*, choosing t'{x,t) : x — e'^^*"* •* we get 
w(t,a;) =e=(*-*')Lj(i',l) ^x"/'' , x > x{t) ^ e"' (23) 

The value uj in this region is therefore determined by the 
boundary; it is a power-law function of x and does not 
depend on time. 

For smaller x, the choice t' — gives 

uj{t, x) = e^^LUQ {xe""^) , x < x{t) (24) 

The influence of the boundary has not spread to this 
inner region yet, and the profile of uj is still determined 




FIG. 1: Illustration to Eq.-s (23), (24): dependence of vorticity 
distribution on time. 



by initial conditions. So, there is no singularity after 
finite time, despite the presence of a power law (Fig.l). 
As time passes, the size x of the inner region decreases, 
and the vorticity profile approaches the singularity but 
never reaches it. 

What happens if there is another boundary condition? 
Substituting arbitrary boundary condition a;(t, 1) = f{t) 
and t'{t, x) : X = e'^^*^* \ we have 



Lo{t, x) = x'/'^f ( t - - Inx ) -^t^^ x^^^it) 



for any given x. Thus, any reasonable (i.e., slower than 
exponential) function / does not change the power law 
and affects only the coefficient, which becomes time- 
dependent (see Fig. 2). Since the boundary conditions 
correspond to large scales, the characteristic time for f{t) 
is of the order of the largest eddy turnover time. This 
statement is also valid for the general case (18). 



Evolution of spectrum 



The solution (23), (24) is not stationary: there is always 
a narrowing non-stationary region X < x{t). We now con- 
sider this solution in terms of Fourier transform. This is 
useful to understand the spectrum evolution (which gives 
the basis to the idea of cascade), and to take account of 
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FIG. 2: Evolution of vorticity distribution in one particular 
case: a = -3, 6 = 1, ujo{x) = 1/ (^1 + [a; + 0.1 sin(107ra::)]^'^^^ , 

uj{t,l) = 1/ + [1 + 0.1e-='*sin(107re^*)]^''^^. One can 

see that the range of strong oscillations drifts to smaller x, 
while inside the 'inertial range' the fluctuations become neg- 
ligible and the power law dominates. 

viscosity. The Fourier transform of vorticity is:^ 

x{t) oo 
x{t) 

The first integral can be rewritten as 



It depends weakly on k for all fc < fc = e~°*, and de- 
creases exponentially as a function of time. The second 
integral is 

oo 
kx{t) 

It is a power law for k < e^"* and decreases sharply at 
larger k. 

So, k) is a step function of k, with the step running 
exponentially to the right as time goes. 

To illustrate this, we consider one particular solution 
of (20), for which the Fourier transform is easy to count 



[7] In fact, (20) is defined for x < 1. But we are interested in 
fc 3> 1 , so there is no difli'erence between Fourier series and Fourier 
transform. 



analytically. This time we do not demand the 'strong' 
boundary condition (21); as we have found in previous 
subsection, it is enough for uj{t, 1) not to grow exponen- 
tially. Let initial distribution of vorticity be 

uJo{x) = {l + ixf''' + {l~ixf''' 

In accordance with (20), (22) evolution of ijj{t,x) takes 
the form: 



L0{x, t) 



1 



■ —at 

le x) 



= 2e^*(l-f a;2e-2-*)^/''^cos (0c/a) 
where tan — xe^°*. For x ^ e"*, we have (f) ~ 7r/2, 

The Fourier transform of this function is: 



uj{k,t) = |fc| 



f'/ag-|*:|e° 



The spectrum falls exponentially at k • 



(25) 

1 = e^"*. The 

result is similar to the effect of viscosity, but the cutoff 
moves along the k axis towards larger values of k (in case 
of dissipation the cutoff would not depend on time) . 

Such a step spectrum spreading to larger k is usually 
interpreted as a cascade, or breaking of vortices. We 
see that in our approach it appears without a cascade, 
and energy is transported to smaller scales by means of 
the narrowing transition region near one selected point, 
which is to become a singular point at infinite time. 

It is usually assumed that viscosity is necessary to get 
stationary statistical picture in turbulence. Indeed, one 
needs viscosity to make all statistical averages, e.g., struc- 
ture functions of all orders, stationary: energy injected 
into a flow at large scales has to be dissipated at viscous 
scale. 

However, our example shows that, in some cases, sta- 
tionary spatial probabilistic distribution can be reached 
even without dissipation in some finite range of scales. 



Effect of viscosity 

It is easy to generalize (19) to include viscosity. Since 
Vu is directed along x axis, the equation takes the form: 

du{x,t) du{x,t) , , . d^u 

ax — h bu(x, t) = J^TT^ 

OX'' 



dt 



dx 



Similarly, the viscous term should be added into the 
right-hand side of (20). Changing to the new variable 

q = xe~°'*, we get 

duj(q,t) , . -yntd^^ 
^-M,,.) = .e— ^ 

(We recall that a < 0.) The Fourier transformation gives: 
a;(A:,t) = e-''*^o(A:e''*)e^'='(i-«"') 
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It appears that, while non-stationarity produces a 'step' 
(exponential fall) running to the right with exponen- 



tial speed k 



viscosity produces a similar (but 



sharper) step which runs to the left, and very quickly 
(after t ~ l/2|a|) becomes stationary at A: ~ \j2\a\lv. 

For the particular example of initial condition consid- 
ered in the previous subsection, instead of (25) we then 
get 



a;(fc,i) = |fc|*/"e-l'=l^"*e^ 



We see that the 'non- viscous' solution does not differ from 
the 'viscous' solution in the range k < e~°*, k < y/2\a\/iy. 

Although the non-viscous equation does not 
'smoothen' the initial perturbations, it transports 
them to smaller and smaller scales and multiplies by 
a decreasing term. So, the solutions of Euler and 
Navier-Stokes equations behave similarly in the limit 
t ^ X. V ^ 0. In this sense, the Euler equation can be 
treated as the inviscid NSE. 



INTRODUCTION OF STOCHASTICS 

In the previous section, the large-scale velocity fluctua- 
tions are treated as deterministic ones. Wc; have seen that 
this causes a power-law dependence of vorticity. This 
would provide a scaling dependence of velocity structure 
functions, but it would be mono-fractal, instead of mul- 
tifractal: the scaling exponents would be proportional to 
their numbers. The multifractal picture can be restored 
if wc take stochastic behavior of large-scale fluctuations 
into account. 

The Theorems cited in Section 3 claim that the 'sys- 
tematic' part (13) of randomly changing matrix (11) not 
only has exponentially growing averages, but also fluctu- 
ations of these exponents are Gaussian random processes. 
An accurate analysis of the stochastic equation (7) with 
account of these fluctuations could probably allow to con- 
struct a complete theory. In this paper, however, we re- 
strict ourselves by a simplifled example from the previous 
Section. 

According to the Theorems, the stochastic generaliza- 
tion of (20) has the form: 

^ + {a + i,{t))x^-{c + Ut))^ = ^ (26) 

Here ^i(i) and £,2{t) are Gaussian random processes with 
zero averages. 

All the relations of the previous Section can be rewrit- 
ten for this case; e.g., (22) becomes 



c{t-t')+J ^2(t")dt" I -a{t~t')-f^i{t")dt" 

co{t,x) = e *' CO \ t',xe 



For a; = 0, taking t' = 0, we get 

Ct+f i2{t")dt" 

w(i,0) = e ° a;(0,0) 

Let ^1 , ^2 be delta-correlated with dispersions Di and 
D2. (For more general case of not delta-correlated pro- 
cesses see Appendix.) The probability density can then 
be written as 

rfP[Ci(t),6(i)] = e lld^i{m2it) 



Thus, 



(w(i,0)") = e""'(e )w"(0,0) 



1 



= 6"^=*+" '°^*/2^"(0,0) (27) 

We see that the averages diverge exponentially as a func- 
tion of time. This characterizes the solution inside the 
non-stationary inner region (24) with growing vorticity. 
The width x of the non-stationary region is determined 
by the condition 



xe 



-at— f ^idt ^ 



But, since / ^idt oc \/i after long time, we have at ^ 
J ^idt and x ~ e"*. 

Thus, adding stochastic fluctuations to a and c in- 
creases the central growth of vorticity but does not 
change significantly the size of the non-stationary region. 

To understand the statistical properties, however, we 
are interested in the outer region x > x{t). In this case, 
by analogy with (23), we choose t'{x,t) in such a way 
that 

t 

Inx = a{t -t')+ j ^idt" (28) 
t' 

Then 

c{t-t')+} i2dt" 

Lj{t,x)=e *' io{t',l) (29) 

The value t'{x,t) is now a random process. An accurate 
calculation of averages of (29) is very complicated, so we 
just make an estimate. We restrict ourselves by small x, 
so I In a; I » 1. An average of / $,\dt" is zero, so we can 
estimate it by / £,\dt" oc ^/t — t' . Thus, in (28) this term 
is much smaller than a{t—t') and can be neglected. From 
(29) we then have 

, / i2dt" 

uj{t,x)c±x''/''e'-'' uj{t',l) 
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-2a 



-2a 



U(x) 




® © o 

^0 ® © (7) 



U(x) 




U(x) 




a-l=(3=-fa/o 



a 



FIG. 3: Cascade model (a) vs Infinite-time singularity (6): tlie same spectrum is produced by different pliysical processes 



Raising to the power n and taking an average, we get 



(UJ ) = X 



dt" 



c 2 Do i/\ c I 2 



nrf6(tV"(i',l) 
(30) 



This scahng of vorticity moments is equivalent to velocity 
structure functions with nonlinear scaling exponents: 

(A«"(/)) ^ (^")/" ~ Z?" , Cn = --n + (31) 

a 2a 

The obtained relations provide an explanation of the non- 
linear dependence of scaling exponents on their order. 



DISCUSSION 

In this paper we derive existence and properties of 
vortex filaments (high-vorticity regions) on the basis of 
the stochastic Eq.(7) (which in turn is derived from the 



NSE) by using the Theorems (page 4). The main result 
is the scaling behavior of vorticity (18), (23) inside these 
vortex filaments, without suggesting finite-time singular- 
ities, and multifractal behavior of vorticity (30) and ve- 
locity (31) statistical moments. Thus, in our approach 
the direct, not only probabilistic formulation of the Mul- 
tifractal model is valid, however there are no singulari- 
ties in the fiow: at any finite time peaks of vorticity are 
smoothed inside constantly narrowing non-stationary re- 
gions. Unfortunately, the numerical values of the coeffi- 
cients are not defined by the theorem. They depend on 
the properties of the random large-scale fiuctuations, and 
a special investigation is needed to derive them. 

The assumptions and simplifications used throughout 
the work are rather general and do not seem crucial. We 
restricted our consideration by symmetric large-scale ve- 
locity gradients Aij, supposing that the large-scale vor- 
ticity can be neglected inside the vortex filament: the 
solution can probably be generalized for all Aij. The 
gaussianity of large-scale fiuctuations is not a necessary 
condition: it is not required in the Theorems, while the 
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Theorems state that fluctuations of the exponents (13) 
are Gaussian. 

Even the simphfled 'straightened' model, apart from 
its illustrative functions, can be valid in the high-vorticity 
regions: the rotation matrix s in (12), being a large-scale 
value, has a characteristic time of changes Tcor '> UJ~^, 
and hence its rotation can be treated as an adiabatic 
process. 

One principal point is separating smaller from larger 

scales, which is not needed in the formal introduction 
of probabilistic description but is used in the 'physical' 
approach (Section 2). The possibility of such separation 
is proved by the mechanism that provides the scaling. 
In the 'canonical' cascade interpretation of the power- 
law spectrum, vortices of all scales are presented and 
contribute to the resulting scaling (Fig. 3a). In this 
approach it is difficult to separate the scales. To the op- 
posite, our consideration shows that the same spectrum 
is produced by small regions near some 'almost singu- 
lar' points (Fig. 3b). Excluding these regions would cut 
the spectrum up to k ~ and the scales are easy 
to separate. An evidence for this second approach comes 
from the numerical simulations [12]: the observed 'coher- 
ent structures' with high vorticity, i.e. vortex filaments, 
are found to be very stable, the lifetime exceeding many 
times the largest-eddy turnover time. This contradicts to 
the idea of cascade. In [12] it is shown that these small- 
scale structures are responsible for the 5/3 law. Picking 
them out breaks the power-law energy spectrum in the 
whole inertial range. 

Returning to the higher-order structure functions, the 
relation (31) proves their power-law dependence and pro- 
vides an explanation of the nonlinear dependence of scal- 
ing exponents on their order. As it was shown in [14, 15], 
the quadratic nonlinearity describes very well velocity 
scaling exponents observed in experiments and numeri- 
cal simulations. More accurate analysis of the stochastic 
equations, with account of rare events which are of most 
importance for high-order structure functions, would of 
course add higher degrees to the expression. But even 
this simplified consideration appears to be enough to 
show that average large-scale exponents Xi in (13) de- 
termine the scaling ('fractal') behavior of the solutions, 
while fluctuations of these exponents produce 'multifrac- 
tality'. 

CONCLUSION 

Thus, in the paper we study the solutions of the 
Navier-Stokes equation in the regions of high vorticity 
(vortex fliaments), treating the large-scale velocity fluc- 
tuations as independent stationary random process. The 
stochastic equations (4), (7) are thus the main equations 
of the paper. 

We analyze the long-time asymptote of the solutions 



to Eq. (7) and show that an inflnite-time singularity 
appears in the limit — > 0, / — > 0; for any finite t, there 
is no singularity, and for any finite I inside the inertial 
range, the solution becomes a power law (18) after some 
time t{l). 

We show that the solution corresponds to random ro- 
tation and systematic exponential stretching of a vortex 
filament. This exponential stretching causes the power- 
law distribution of vorticity, the resulting spectrum is 
quite similar to that expected from the model of break- 
ing vortices. 

Taking into account the stochastic component of the 
stretching, wc get the multi-scaling distribution of vortic- 
ity (and velocity differences), and quadratic dependence 
of velocity scaling exponents on their order (31). As it 
was shown in [14], this result agrees very well with ex- 
perimental and DNS data. All the results do not depend 
on the assumptions on the properties of the large-scale 
random process. 

We thank Prof. A.V. Gurevich for his permanent in- 
terest to our work. The work was partially supported by 
RAS Program 18 IL. 

APPENDIX: CALCULATION OF STATISTICAL 
MOMENTS IN THE CASE OF FINITE 
TIME-CORRELATED RANDOM PROCESS 

The averages (27), (30) are calculated under the as- 
sumption of delta-correlated coefficients ^1,^2 in Eq. 
(26). However, the same result can be obtained in the 

limit t ^ 00 without the assumption. To illustrate this, 

t 

k f a{ti)dti 

we consider the average (e ° ) where a{t) is a ran- 

dom Gaussian process with correlation function 

00 

{a{h)a{t2)) = G{h -t2) , J dTG{T) = 1 (32) 

— 00 

The probability density of a can be written in the form 
[22]: 

t t 

f dti f dt2G-'^{ti-t2)a{ti)a{t2) 

P[a](t) = e ° 
where is defined by 

/oo 
dt'G{t2 - t')G-\t' - h) = 6{t2 - h) (33) 
-CSO 

For the statistical moments we then get 

t 

k f a(ti)dti 

(e ) = (34) 

t t t 

f -5/^*1 J dt^G-^ (H-t2)a(ti)a(t2)+k J dtia(ti) 

= / M da(T)e 00 
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This is a Gaussian integral, thus, the saddle-point 
method gives an exact result in the case [22]. The opti- 
mal trajectory is defined by the condition: 



(t t 
-i I dh [ dt2G-\h - t2)a{ti)a{t2) 



+k 







J dtia{ti) 



Hence, 



j dtiG~'^{t' -ti)ao{ti) = k , 0<t'<t (35) 



Making use of (33), we get 



oo(i') = k j dtiGii! - h) 



(36) 





Substituting to (34) we get 

i 



k J a{ti)dti 

(e ° ) ~ exp 



ao{ti)dti X 



-i J dt2G-\h - t2)ao{t2) + k 



and with account of (35), (36): 



kf a(ti)dti 

{e ° ) ~ exp 



= exp 



k / dtiao{ti) 



(37) 



^k^ j dti j dt2G{ti - 12) 







(One could as well obtain the same result by shifting a{t) 
to get the perfect square in the exponent.) 

Prom (37) it follows that as f ^ oo (or, more precisely, 
for any t larger than the correlation time) the statistical 



moments grow exponentially, just as 
correlated process. 



in the case of delta- 
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